Upper limits on gravitational waves emission in association with the Dec 27 2004 

giant flare of SGR1806-20 
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At the time when the giant flare of SGR1806-20 occurred, the AURIGA "bar" gw detector was 
on the air with a noise performance close to stationary gaussian. This allows to set relevant upper 
limits, at a number of frequencies in the vicinities of 900 Hz, on the amplitude of the damped gw 
wave trains, which, according to current models, could have been emitted, due to the excitation of 
normal modes of the star associated with the peak in X-rays luminosity. 



PACS : 04.80.Nn, 95.55.Ym 

On 27 December 2004 the Soft Gamma-ray Repeater 
SGR1806-20 gave a giant flare, which was observed by a 
number of instruments Q . 

The flucnce, if the emission is assumed isotropic, at the 
distance of d ~ 15 kpc would imply an energy some hun- 
dred times larger than any other known giant flare 0, 13 • 
Soft gamma-ray repeaters are thought to be magnctars 
(see and refs. therin). It has been suggested that 
the extreme energy event of 27 December 2004 is due to 
a catastrophic instability involving global crustal failure 
and magnetic reconnection Observations by CLUS- 
TER and TC-2, in combination of data from GEOTAIL, 
gave evidence that the steep initial rise contains two ex- 
ponential phases, of e-folding times 4.9 ms and 67 ms 
respectively, which covered the 24 ms before the time 
of the peak intensity t p ; all the timescales support the 
notion of a sudden reconfiguration of the stars magnetic 
field, producing large fractures in the crust jfj . In partic- 
ular these authors remark that the intermediate sa 5 ms 
time is naturally explained if the rising time is limited by 
ithe propagation of a triggering fracture of size w 5 km, 
as it would be predicted by the theory of reference [|J . 

According to a few somewhat different models, as a 
consequence of crustal cracking Q or reconfiguration of 
the moment of inertia tensor |8j , non-radial kHz oscilla- 
tion modes of the neutron star would be excited, giving 
emission of gravitational waves (gw), possibly at frequen- 
cies where the gw bar detector AURIGA jf| is sensitive 



(see insert of Fig. 1). Both the above quoted mod- 
els predict gw emission, starting very close to t p , which 
involves kHz iron radial modes of oscillation of a neu- 
tron star with few hundred ms damping time. The ex- 
pected waveforms can be approximately parametrized as 
h(t) = ho exp(— t/Tg) sin(27r/ s t), where ho is the maxi- 
mum gw amplitude, f s and t s are the frequencies and 
damping times of normal modes; the polarization of the 
wave is not known. The frequencies of the various modes 
are still under study and depend on a variety of factors 
as EOS, temperature, density, age, rotational state of the 
star, etc. |lf| so that we are unable to anticipate with 
any confidence what specific set of gw emission frequen- 
cies could be the one expected for a magnetar ready to 
undergo a supcrgiant flare. Still the lowest lying modes, 
g— , /— and, marginally, p — modes, could well be in the 
frequency range 500 -f- 1500 Hz, depending on the status 
of the star. 

Within a factor 10 in gw amplitude, AURIGA is sen- 
sitive to gws from ~ 800 to 1050 Hz. Here we limit the 
analysis to the most sensitive part of the band, namely 
between 850 and 950 Hz (see the insert in Fig. 1), where 
the detector sensitivity varies no more than a factor 4 in 
amplitude. Since the last upgrading of the suspensions 
on Dec 2 nd 2004 the detector is well behaved in the sense 
that performs stationary gaussian, after epochs of envi- 
ronmental disturbances are vetoed by means of auxiliary 
channels (i.e. signals at frequencies where the detector is 
gw insensitive). During nights and week-ends the vetoed 
epochs becomes less frequent and shorter, so that the 
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detector achieve close to 90% stationary gaussian opera- 
tion; in particular on time spans of minutes we can use 
the data, without even applying vetoes. This is the case 
for the time span of about ±100s around the epoch of 
the 27 December 2004 giant flare of SGR1806-20, which 
we use in this analysis. We show in the following that the 
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FIG. 1: Plot of £ 1/2 in the frequency band 930 935 Hz as 
a function of time; the origin in the x-axis corresponds to the 
arrival time of the flare X of SGR1806-20 at the AURIGA site. 
The insert shows the AURIGA one-sided noise spectral den- 
sity, as expressed in equivalent gw amplitude h r (t) at input. 
The vertical dashed area shows the position of the frequency 
bin 930 Hz (see text). 

noise is driven by a zero mean stochastic gaussian process 
with a stationary correlation function. For what concern 
the directional sensitivity, the orientation of AURIGA in 
respect to the direction of SGR1 806-20 was such that 
the antenna pattern, averaged over polarizations, gave 
maximal sensitivity at the time of the giant flare. Then 
we have a unique opportunity to search in our data for 
gravitational waves emitted at the peak time of the gi- 
ant flare. We take the peak time t p to be 21:30:26.68 
UT of 27 December 2004 after taking into account the 
time difference between the arrival time at INTEGRAL 
[HI and at AURIGA site of 133.427 ms.0 This time 
corresponds also to the peak position of the CLUSTER 
data which show, after the last exponential rise, the evi- 
dent start of a phase in which damping occurs until the 
signal gets below 1/10 of the peak value, ~ 300 ms after 
the peak. Q Following the models quoted above, in both 
cases we can assume the peak time t p as the start of the 
gw excitation and t s = 100 ms, that is 1/3 of 300 ms, 
as the corresponding damping time. In order to extract 
the signal power first we reconstruct the gw amplitude 
h r {t) at input through the detector transfer function. 
Then we slice the gw sensitive frequency band of AU- 



RIGA in contiguous and non-overlapping sub-bands fj 
of constant width Af, and centered in /J = fj + A//2, 
by means of digital top-hat filters in the frequency do- 
main Tj(f) = #(|/| — fj) — #(|/| — fj — A/). Within each 
sub-band, we compute the the equivalent input signal 
power over a time span At 



£, 



Tj * h 2 r {t + kAt) dt 



(1) 



A/ 



where * stands for time convolution. The £j(t) is sam- 
pled every At to construct the time series £j(k) with 
k integer. We decided a priori a fixed partition of 
the time frequency plane: A/ = l/(2r s ) = 5 Hz and 
At = 201.5 ms ~ 2t s . For each sub-band fj, we ana- 
lyzed the resulting time series of £j (k) over a time span of 
±100 s around the peak time t p to check the "off source" 
noise statistics. The £j{k) sample including the peak 
time t p is then compared to the measured noise statis- 
tics, looking for any evidence of excess power. To be more 
precise, the a priori choice of our sampling time made t p 
to fall 120 ms after the beginning of the integration time 
At of the "on source" sample. Fig. 1 shows how £j fluc- 
tuates on the time spans of ±5 s around the time of the 
flare t p for the sub-band fj = 930 Hz. A gw emission 
at frequency f s would give an excess power in the band 
A/ centered at the // such that \f s - f?\ < Af/2. The 
released energy would be maximum in the "on source" 
sample. The excess signal power in each sub-band Af 
can be easily calculated from the expected waveform and 
reads 
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where x = (2 7 rr s )- 1 /A/ and 6 = \f s - f?\/Af are the 
ratios between the signal bandwidth (= l/(27rr s )) and 
the detuning of the signal frequency and the bandwidth 
Af, respectively. With our choice of parameters for the 
analysis, the excess signal power is approximately (within 
a few % error) 
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where A/ e // = 4 Hz. To check the statistics of 
the "off source" samples, we histogram each time se- 
ries £j(k) and compare them with the predicted prob- 
ability density functions assuming gaussian noise, by 
fitting for the variance separately in each sub-band. 
The fitting probability density function is a x 2 distri- 
bution with a effective degrees of freedom p(£;a 2 ) = 
2- Q / 2 (f/a 2 ) Q / 2 - 1 exp(-£/2a 2 )/r(a/2)/cT 2 , where a 2 is 
the variance of the underlying gaussian stochastic pro- 
cess. We show in Fig. 2 the close agreement with predic- 
tion of the data for the frequency bin fj = 930 Hz, over 
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FIG. 2: Histogram of £ in the frequency band 930 -j- 935 Hz; 
solid line represents the best fit curve of a x 2 distribution 
with a = 3.6 effective degrees of freedom. The insert shows 
the p-level distribution of the same fit for the histograms of 
the 18 frequency bins. 



the extended time span of ±100 s. The results for all 
other frequency bins are similar. The goodness of the fit 
has been checked by a x 2 test, and the resulting p-values 
for all the sub-bands are consistent with a uniform dis- 
tribution in the unit interval, as expected (see the insert 
of Fig. 2 and Tab.l). In Table 1 we report the parameter 
a 2 of the fit of p(£) to the experimental data. We find 
that the dependence on the effective degrees of freedom 
a of the p-levels is very weak and, within the statistical 
errors, we can fix a = 3.6 for all the frequency bins. The 
p-level distribution is uniform in the unit interval (see the 
insert of Fig. 2). The stationary behavior, at least for 
timescales of few minutes, is shown by the constancy in 
time of the parameters needed to fit the noise model. We 
take advantage of the classical theory of hypothesis test- 
ing to establish if the samples £ tp corresponding to the 
arrival time t p are affected by the presence of a gw signal. 
To test the null hypothesis Hq, i.e. that the sample is 
drawn from the estimated noise probability distribution 
in absence of signals, we set a threshold £ cr correspond- 
ing to a confidence level (C.L.) p{£ < £ cr ) > 1 — p cr - 
The threshold for 1 — p cr = 95 % C.L. corresponds to 
£ cr = 8.8 x a 2 . Thus one sees from Table 1 that no excess 
of gw power is found at t p and therefore we have to set 
up upper limits. We set conservative confidence intervals 
for £ s using a confidence belt construction which en- 
sures non-uniform coverage greater or equal to 90%. The 
confidence belt construction proceeds as follows. Assume 
that the signal magnitude is £ s . The measured £ in each 
sub-band (Eq. 1) obeys a non-central % 2 distribution 
with central parameter equal to £ s /(t 2 (here we drop for 
simplicity the index of the sub-band). Its corresponding 



TABLE I: List of fit parameter a of the histogrammed £ data 
samples ±100 s around t p tabulated as increasing sub-band 
frequencies; the data for the 870 Hz sub-band have been dis- 
carded a priori as this band is contaminated by environmental 
noise. The "on source" value of £ including the trigger time 
t p is also reported as well as the computed upper limit with 
confidence > 95%. 
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probability density function can be written as: 



p(£;£ s ,a) = ^ exp ( - 



£ + £ s \ (£_ 
2a 2 ) U 



(a-2)/4 



4/2-1 \j££slo 2 ) , 



(4) 



where Ik (x) are the modified Bcsscl functions of the first 
kind of order fc. The g-quantile of this distribution, 
£ q (£ s ,a), is implicitly defined by q = J Q q p(£;£ s ,a)d£. 
For each value of the unknown £ s we define the 95% con- 
fidence belt boundaries £ hi and £i ow as 



£hi{£ Sl o) 



if £ s < £f{cr) 
£^%(£ s ,a) otherwise (5) 

£ low {£ s ,&) = £95% {£s,<?) 



where £f is implicitely defined by £^%{£f , er) = 
f 95 % (0, a). This confidence belt defines a set of confi- 
dence intervals on £ s , whose frcquentist coverage is - by 
construction - 90% for £ s > £f, and 95% for £ s < £f. 
In other words, for every value oi£t p from each sub-band, 
if £t p < £95% (0, Oj) we set an upper limit equal to £f , 
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otherwise our procedure gives a two-sided confidence in- 
terval. In all sub-bands we obtain upper limits, which can 
be written as £f ~ 18 X er|. These limits ranges from 

S 1 / 2 = 3.5 x 10-21^-1/2 t0 £ 1/2 = L4 x IQ-^Hz' 1 ' 2 , 

according to AURIGA sensitivity. 

The initial amplitude of the neutron star normal modes 
ho is related to £ by Eq. (3) that gives, for the best upper 
limit, h < 2.7 x 1CT 20 . We discuss now the upper limit 
in terms of the total gw energy e gw = E 9W /Mqc 2 emitted 
by the normal modes excitation during the peak of the 
giant flare of SGR1806-20. The well known formula of the 
quadrupolar radiation, for the expected gw signal 01 , can 
be written as ho = (e gw ci?s/(47r 2 r s ) 1 / 2 / (f s d), where Rs 
is the Swartzchild radius of one solar mass black hole. 
Thus the resultant upper limit on e gw reads 



c-gw 



< 3 x 10" 



1.3 x lO- 4 ^- 1 

2 



15 kpc\ /930 Hz 



We should notice that a gw bar detector has a polariza- 
tion dependent sensitivity; hence, for an unpolarized or 
linearly polarized gw, the result in Eq. (6) should be 
multiplied by a factor 2 or cos 2 (2tp) respectively, where 
?p is the angle between the bar axis and the polarization 
of the wave. We conclude that, if the star ever emitted 
gws from excitation of its normal modes at any of the 
frequencies studied here, in the time span At containing 
the flare time t p , the gw amplitudes and energetics are 
limited as above. If the giant flare of SGR1806-20 on 27 
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FIG. 3: Upper limts on the gw energy released at the source 
around the flare peak time t p expressed as a fraction of Mqc 2 . 



December 2004 is indeed some 100 times more energetic 
(however see ref. [l^l and if the gw luminosity scales 
with the em luminosity, then, for the frequencies consid- 
ered, our upper limits come close to the predictions of 
the models of refs. 0, @ which give an energetics of the 
order of e gw ~ 5 x 10 -6 . The method used here is of 
course sub-optimal and the upper limits are somewhat 
weaker than the "optimal" matched filter. In any case 
this work shows that, as there is the specific peak time 
t p to be used as external trigger, it is worth to make 
searches even with a single detector if its noise is well 
behaved. An extension of such searches involving the gw 
detectors on the air in a coincidence search, would also 
allow to use the information of the gw travel delays be- 
tween the detectors to select against spuria, and would 
give the most exhaustive and efficient search, in terms 
of frequency coverage and confidence in improving the 
limits, if not to get a candidate detection. 



We are gratefull to Roberto Turolla for a critical read- 
ing of the manuscript. 
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